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Abstract 

First-principles  density  functional  theory  (DFT) 
calculations  using  the  generalized  gradient 
approximation  (GGA)  have  been  conducted  to  study  the 
adsorption  of  a series  of  high-nitrogen  compounds  of 
increasing  sizes  and  complexity  on  the  Al(lll)  surface. 
The  calculations  employ  periodic  slab  models  with  4 Al 
layers,  ranging  in  size  from  (3*3)  to  (7*7)  surface  unit 
cells,  and  containing  up  to  196  Al  atoms. 
Complementary  quantum  chemical  calculations,  utilizing 
DFT  and  second-order  perturbation  theory  methods,  of 
the  ground  state  potential  energy  surfaces  of  the 
corresponding  polynitrogen/high  nitrogen  species  in  the 
absence  of  the  aluminum  surface  also  have  been 
performed.  For  the  set  of  chemical  species  Nx(x=l,5), 
NHx(x=l-3),  N2Hx(x=1-4)  and  N3H,  N3H3,  and  N4H4,  the 
adsorption  configurations  at  different  surface  sites  and 
the  corresponding  binding  energies  have  been 
determined.  This  analysis  has  been  further  extended  to 
high-nitrogen  compounds  N5H  and  NAT?.  For  these  two 
systems  it  was  found  that  the  initial  bonding  to  the 
surface  takes  place  through  a molecular  mechanism 
(nondissociatively)  with  involvement  of  single  or  multiple 
N atoms  of  the  molecule.  However,  dissociation  on  the 
surface  can  take  place  with  small  activation  energies. 
This  set  of  calculations  has  been  further  extended  to 
include  1,3,5-triazene  and  1,2,3-triazine  (C3H3N3)  as  well 
as  larger  substituted  triazene  systems  such  as  C9N30  and 
CisNls.  For  these  large  systems  it  was  found  that 
bonding  takes  place  through  multiple  N centers  with 
formation  of  highly  strained  and  deformed  adsorption 
configurations.  In  a number  of  instances  the  adsorption 
takes  place  dissociatively  with  N2  elimination. 
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1.  Introduction 

Investigation  of  the  interaction  of  energetic  high- 
nitrogen  compounds  with  ultrafine  or  nanophase  particles 
of  aluminum  are  a topic  of  current  interest  with  potential 
Department  of  Defense  (DoD)  applications  in  rocket  and 
missile  propulsion.  One  of  the  topics  of  prime 
importance  is  to  determine  if  the  thermal  instability 
inherent  in  many  high  nitrogen  compounds  might  be 
reduced  by  adsorption  or  chemisorption  onto  the  surface 
of  aluminum,  while  simultaneously  forming  a protective 
coating  on  the  metal  surface  which  minimizes  the 
formation  of  an  inert  oxide  surface  layer. 

2.  Computational  Method 

The  calculations  performed  in  this  study  were  done 
using  the  Vienna  ab  initio  simulation  package 
(VASP)11'3'.  This  program  evaluates  the  total  energy  of 
periodically  repeating  geometries  based  on  density- 
functional  theory  and  the  pseudopotential  approximation. 
In  this  case  the  electron-ion  interaction  is  described  by 
fully  non-local  optimized  ultrasoft  pseudopotentials 
similar  to  those  introduced  by  Vanderbilt14,51.  Periodic 
boundary  conditions  are  used,  with  the  one-electron 
pseudo-orbitals  expanded  over  a plane-wave  basis  set 
with  a cutoff  energy  of  435  eV.  Calculations  were 
performed  using  the  generalized  gradient  approximation 
(GGA)  density  functional  theory  with  PW91  exchange- 
correlation  functional161.  The  sampling  of  the  Brillouin 
zone  was  performed  using  a Monkhorst-Pack  scheme171. 
The  minimum  energy  path  between  different  minima  was 
determined  using  the  climbing  image  nudged  elastic  band 
method  developed  by  Jonsson  and  co-workers18,91. 

The  GAMESS[10]  quantum  chemistry  code  was  used 
to  predict  the  structures  and  properties  of  isolated 
polynitrogen  and  high-nitrogen  compounds.  The 
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theoretical  methods  employed  included  second  order 
perturbation  theory^1 11  (MP2,  also  known  as  MBPT(2)) 
and  density  functional  theory  (DFT)  using  the  hybrid 
B3lyP[12]  functional.  The  6-31  lG(2df,p)  basis  set[13] 
was  used  throughout. 

3.  Results  and  Discussion 

In  order  to  provide  a comprehensive  analysis  of  the 
interaction  of  polynitrogen  systems  with  aluminum 
surface  in  our  studies,  we  focused  on  four  different 
categories  of  high  nitrogen  and  polynitrogen  systems  of 
increasing  complexity.  The  first  set  is  composed  by 
individual  atomic  and  molecular  all  nitrogen  systems 
from  the  series  Nx  (x=l,5).  The  second  set  is  formed  by 
small  hydrogenated  nitrogen  systems  including 
NHx(x=l,3),  N2Hx(x=1-4),  N3Hx(x=1,3),  N4H4,  N5H,  and 
N6H2.  The  third  set  includes  systems  of  increased 
complexity  such  as  (1,3,5)  and  (1,2,3)  triazene  (C3N3H3) 
as  well  as  mono-substituted  derivatives  of  triazenes  such 
as  C7N8H2  and  C5N,2H2.  Finally,  the  fourth  set  is 
represented  by  large  high-nitrogen  compounds  of 
trisubstituted  triazene  molecules  such  as  C9N30  and 
CijN18. 

Given  the  large  variety  of  molecular  sizes  involved, 
the  adsorption  studies  of  the  above  compounds  were  done 
using  slab  models  varying  from  (3><3)  for  the  case  of 
small  molecular  systems  to  (7x7)  surface  units  for  the 
largest  systems.  In  all  cases  the  slab  models  contained 
four  layers  of  atoms.  The  optimizations  have  been 
performed  for  both  the  adsorbate  and  the  A1  atoms  in  the 
top  two  layers  while  the  bottom  two  layers  were  frozen  at 
the  bulk  optimized  conditions.  Adsorption  at  different 
surface  sites  and  for  different  molecular  orientations  have 
been  determined  and  the  corresponding  adsorption 
energies  have  been  calculate  using  the  expression 
Eads- Em0[ec  + Esiab  — E(:iKm:C- slab),  where  Em0]ec  is  the  energy 
of  the  isolated  molecular  species  at  its  equilibrium 
geometry,  ES[ab  is  the  total  energy  of  the  isolated 
aluminum  slab,  and  E(moiec+siab)  is  the  total  energy  of  the 
adsorbate/slab  system.  A positive  Eads  corresponds  to  a 
stable  adsorbate/slab  system.  In  the  present  study,  only 
neutral  species  have  been  considered. 

Nj  (x=l,5)  Adsorption  on  Al(lll).  The  binding  of 
individual  Nx(x=l,5)  has  been  determined  using  the  (3x3) 
supercell  model  and  representative  adsorption 
configurations  are  represented  in  Figure  1.  The  most 
stable  species  was  found  to  be  the  N atom  which  binds  to 
the  surface  in  a three  fold  configuration  with  an  energy  of 
150.2  kcal/mol  at  a face-centered  cubic  (fee)  surface  site 
and  of  142.2  kcal/mol  at  a hexagonal  closest  packed  (hep) 
site.  By  increasing  the  atomic  coverage,  the  binding 
energy  further  increases  to  161.5  kcal/mol  at  full 


coverage  when  all  fee  and  hep  sites  are  occupied.  In  this 
case,  the  resulting  atomic  configuration  resembles  the 
aluminum  nitride  (AIN)  phase.  In  the  case  of  the  N2 
molecule,  the  adsorption  on  the  surface  is  not  stable 
relative  to  isolated  N2  in  the  gas  phase.  The  N3  species 
can  adsorb  by  simultaneous  binding  to  three  A1  surface 
atoms  and  maintaining  a triangular  gas  phase  geometry 
(Eads=62  kcal/mol)  (see  Figure  lc)  or  in  an  open  chain 
geometry  which  involves  bonding  to  four  A1  surface 
atoms  (Eads=75.7  kcal/mol)  (see  Figure  Id).  In  the  latter 
case  we  note  that  the  fourth  A1  atom  involved  in  bonding 
is  significantly  pulled  up  from  the  surface.  In  the  case  of 
the  N4  molecule,  several  geometries  have  been 
considered,  including  tetrahedral,  rectangular,  and  linear 
geometries.  The  tetrahedral  geometry  was  found  to  be 
either  weakly  physisorbed  on  the  surface  or  unstable  with 
respect  to  dissociation  into  two  N2  molecules.  For  the 
other  two  isomsers  of  N4,  stable  adsorption  sites  were 
found  as  indicated  in  Figures  le  and  If,  where  three  or  all 
four  of  the  N atoms  are  involved  in  direct  bonding  to  the 
surface  A1  atoms.  Finally,  for  the  N5  species  our 
calculations  indicate  that  the  most  stable  configuration 
corresponds  to  a vertical  geometry  where  bonding  can 
take  place  either  on-top  site  (see  Figure  lh)  or  in  a bridge- 
like configuration  (see  Figure  Ig).  In  the  latter  case,  the 
binding  energy  is  90  kcal/mol  relative  to  the  isolated  N5 
species  and  slab. 

The  energies  of  different  adsorbed  Nx(x=l,5)  species 
are  presented  in  a common  diagram  in  Figure  li),  in 
which  the  energy  zero  is  taken  to  be  the  sum  of  the 
isolated  A1  slab  and  that  of  three  isolated  N2  molecules  in 
the  gas  phase.  From  this  diagram  it  can  be  easily  seen 
that  the  adsorbed  N atom  is  the  most  stable  species, 
followed  by  the  N3  and  N4  molecular  species.  This 
diagram  also  indicates  that  for  the  adsorbed  Nx  (x  > 1) 
species,  further  dissociation  leading  to  formation  of 
adsorbed  N atoms  is  energetically  favorable. 

NjHy  Adsorption  on  Al(lll).  The  adsorption 
properties  of  the  series  of  compounds  NHx(x=T-3), 
N2Hx(x=1-4),  N3Hx(x=1,3),  N4H4,  N5H,  and  N6H2  have 
been  predicted  using  the  (3x3)  slab  model  described 
above.  Representative  adsorption  configurations  of  these 
systems  are  presented  in  Figure  2.  Among  the  most 
strongly  bound  species  in  this  series  are  NH  (Eaijs=130 
kcal/mol),  NH2  (Eads=72.9  kcal/mol),  NNH2  (Eads=76.3 
kcal/mol),  and  N3H3  (Eads=75.9  kcal/mol).  The  larger 
systems  in  this  series  have  significantly  smaller  binding 
energies,  with  values  of  14.1  kcal/mol  for  N2H4,  8.7 
kcal/mol  for  N5H,  and  9.9  kcal/mol  for  N^. 

Minimum  Energy  Path  for  N5H  and  N6H2 
Dissociation.  Starting  from  the  most  stable  adsorption 
configurations  of  N5H  and  N6H2,  the  minimum  energy 
pathways  for  molecular  dissociation  of  these  species  have 
been  characterized.  Two  representative  pathways  are 
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shown  in  Figures  2q  and  2r,  respectively.  In  the  pathway 
shown  in  panel  q),  one  of  the  N-N  bonds  of  the  N5H  ring 
breaks  with  a small  activation  energy  of  about  4.8 
kcal/mol.  Following  this  step,  a second  N-N  bond  of  the 
remaining  ring  dissociates  without  any  additional  barrier, 
leading  to  formation  of  two  N2  molecules  and  the  NH 
radical  adsorbed  at  a three-fold  hollow  site.  This  reaction 
is  highly  exothermic,  by  -92.5  kcal/mol.  Similar  behavior 
is  observed  in  the  case  of  N6H2  dissociation  (see 
Figure  2r).  In  this  case,  the  NF12  group  dissociates  from 
the  N5  ring  with  a small  activation  energy  of  1.8  kcal/mol. 
In  the  final  configuration,  both  the  NH2  group  and  the  N5 
ring  are  adsorbed  on  the  surface.  The  overall  dissociation 
reaction  is  also  exothermic,  by  -63.7  kcal/mol.  These 
results  suggest  that  dissociation  of  N5Fl  and  N6H2  on  the 
A1  surface  is  very  efficient. 

Adsorption  of  Triazene  and  Mono-substituted 
Triazene.  Representative  adsorption  configurations  of  a 
(1,3,5)  and  (1,2,3)  triazene  molecule  are  shown  in  Figures 
3a-d,  which  have  been  calculated  using  a (6*6)  surface 
slab  model.  For  (1,3,5)  triazene,  both  adsorption  through 
a single  nitrogen  atom  (Figure  3a)  or  through  multiple 
bonds  involving  both  N and  C atoms  (Figure  3b)  can  take 
place,  leading  to  total  binding  energies  of  11.2  and  16.9 
kcal/mol,  respectively.  In  the  case  of  (1,2,3)  triazene,  two 
or  three  of  the  N atoms  of  the  molecule  can  be  involved  in 
bonding  to  surface,  which  increases  the  binding  energies 
to  37.2  kcal/mol  for  the  q4  configuration  depicted  in 
Figure  3c)  and  3 1.4  kcal/mol  for  the  q3  structure  in  Figure 
3d). 

The  adsorption  configurations  of  mono-substituted 
triazenes  C7N8H2  (see  panels  (e-g)  in  Figure  3)  and 
C5N12H2  (panels  (h-j)  in  Figure  3)  have  also  been 
analyzed.  In  these  instances,  the  adsorption  takes  place 
with  significant  deformations  of  the  molecular  structure 
and  with  involvement  of  different  sites  of  the  molecule. 
In  the  case  of  C7N8H2,  adsorption  takes  place  molecularly 
and  the  overall  binding  energy  is  relatively  small  with  a 
maximum  value  of  6.80  kcal/mol  for  the  structure  in 
Figure  3g.  A significant  portion  of  the  stabilization 
gained  from  formation  of  Al-N  or  Al-C  bonds  is  lost  in 
the  deformation  of  the  gas  phase  structure  of  the  C7N8F12 
molecule.  In  the  case  of  the  C5N|2F12  system,  similar 
interactions  with  the  surface  are  observed,  with  formation 
of  multiple  bonds  (for  example  eight  bonds  are  present  in 
the  structure  structure  shown  in  Figure  3i)  and  significant 
structural  deformation  of  C5Ni2H2.  In  addition, 
dissociative  adsorption  via  N2  elimination  can  take  place, 
as  shown  in  Figure  3,  panel  (j).  Note  that  the  strong 
molecular  interactions  between  the  adsorbate  and  the 
metallic  surface  leads  to  significant  upward  shifts  of  some 
of  the  A1  surface  atoms. 

Adsorption  of  Trisubstituted  Triazenes. 

Calculations  in  this  case  have  been  done  using  the  largest 


supercell  models  employed  in  this  study,  containing  (7x7) 
surface  units  and  4 layers  of  aluminum  atoms.  The 
overall  trends  observed  for  mono-substituted  triazene  are 
also  seen  in  the  case  of  trisubstituted  triazenes  C9N30  and 
Ci5N18  (see  Figures  3k-o).  Specifically,  the  adsorption  to 
the  surface  can  take  place  through  formation  of  multiple 
bonds,  primarily  between  A1  and  N.  As  a result, 
significant  deformations  and  structural  changes  take 
place,  both  within  the  adsorbate  molecule  and  the 
uppermost  layer  of  A1  atoms,  several  of  which  are 
displaced  upwards.  Binding  energies  as  high  as  69.4 
kcal/mol  are  calculated  for  C]5Ni8  (see  Figure  3m).  Note 
that  in  the  case  of  the  C]5N18  molecule,  which  contains  six 
N3  groups,  elimination  of  several  N2  molecules  can  take 
place  (panel  (0)  in  Figure  3.) 

4.  Conclusions 

A systematic  study  of  a series  of  polynitrogen  and 
high  nitrogen  species  ranging  from  small  Nx(x=l,5)  and 
N2Hx(x=l-4)  to  complex  C9N30  and  C15Ni8  systems 
interacting  with  the  aluminum  (111)  surface  have  been 
investigated  using  plane-wave  DFT  calculations  in 
conjunction  to  slab  models  as  large  as  four  layers  of  (7X7) 
surface  units.  Our  calculations  indicate  that,  with  the 
exception  of  N2,  the  majority  of  polynitrogen  and  high 
nitrogen  compounds  will  bind  to  the  surface  through 
formation  of  Al-N  bonds,  while  formation  of  Al-C  bonds 
is  observed  in  only  a few  instances.  For  the  ensemble  of 
molecular  systems  considered,  a significant  range  of 
variation  in  binding  energies  has  been  determined 
depending  on  the  number  of  surface  bonds  involved, 
surface  site  and  molecular  orientation.  Additionally,  in 
the  great  majority  of  cases  investigated,  adsorption  was 
found  to  take  place  with  important  molecular 
deformations  and  often  with  significant  displacement  of 
aluminum  atoms  from  the  surface.  For  the  systems 
containing  N3  groups  (C5Ni2H2  and  C9N30),  adsorption 
leads  to  elimination  of  multiple  N2  molecules.  For 
smaller  systems  such  as  N5H  and  N6H2,  small  activation 
energies  were  found  for  dissociation  pathways  leading  to 
formation  of  NH  and  NH2  species  adsorbed  on  the 
surface. 

Future  work  will  extend  this  set  of  investigations  to 
include  the  thermal  effects  upon  surface  adsorption  and 
adsorption  on  nanoclusters  of  different  sizes. 
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Figure  1.  Side  and  top  views  of  the  most  stable  adsorption 
configurations  of  (a)  N,  (b)  N2,  (c),  (d)  N3,  (e),  (f)  N4  and  (g),  (h) 
Ns  on  the  Al(111)  surface.  Panel  i)  indicates  the 
corresponding  potential  energy  diagram  for  the  relative 
stability  of  the  N«(x=1,5)  species  on  the  Al(111)  surface.  The 
sum  of  energies  of  the  isolated  slab  and  three  N2  molecules 
in  the  gas  phase  was  taken  as  the  zero  of  energy.  The  labels 
(a)  and  (ph)  denote  adsorbed  and  physisorbed  species, 
respectively,  while  (ch),  (tr)  and  (rect)  denote  chain, 
triangular,  and  rectangular  configurations,  respectively. 


